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ABSTRACT 

Magnetically aligned velocity anisotropy over varying physical conditions and envi- 
ronments within the Taurus Molecular Cloud is evaluated from analysis of wide field 
spectroscopic imaging of ^^CO and ^'^CO J=l-0 emission. Such anisotropy is a re- 
sult of MHD turbulence in the strong magnetic field regime and provides an indirect 
measure of the role of magnetic fields upon the gas. Velocity anisotropy aligned with 
the local, projected mean magnetic field direction is limited to fields with low sur- 
face brightness ^^CO emission corresponding to regions of low visual extinction and 
presumably, low gas volume density. The more optically thin ^'^CO J=l-0 emission 
shows little evidence for velocity anisotropy. We compare our results with computa- 
tional simulations with varying degrees of magnetic field strength and Alfvenic Mach 
numbers. In the diffuse, molecular envelope of the cloud, a strong magnetic field and 
sub- Alfvenic turbulent motions are inferred. Super- Alfvenic motions are present within 
the high column density filaments of the Taurus cloud. From this trans- Alfvenic fiow, 
we constrain the scaling exponent, k, of the magnetic field density relation {B n") 
to be near zero as expected for ambipolar diffusion or material loading of magnetic 
fiux tubes. 

Key words: MHD - turbulence - ISM: magnetic fields - ISM: molecules - ISM: 
kinematics and dynamics 
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1 INTRODUCTION 

The origin of the supersonic motions observed within molec- 
ular clouds is a fundamental but poorly understood aspect of 
interstellar matter and star formation. While such motions 
are attributed to turbulent magneto-hydrodynamic (MHD) 
flows within the cloud, it is currently uncertain whether ve- 
locities of such flows are less than (sub- Alfvenic) or greater 
than (super- Alfvenic) the Alfven velocity, va, in the cloud. 
This distinction is important in the context of describing the 
formation of molecular clouds and the interstellar processes 
that regulate the production of stars. Turbulent fragmen- 
tation considers supersonic, super-Alfvenic turbulence that 
produces strong shocks at the interface of converging flows. 
The resultant shock-induced density perturbations are seeds 
to gravitationally unstable cores from which individual stars 
and clusters condense (Padoan & Nordlund 2002; Krumholz 
& McKee 2005). For sub- Alfvenic clouds, the observed mo- 
tions are attributed to a spectrum of Alfven waves (Arons 
& Max 1975; Kudoh & Basu 2003; Mouschovias, Ciolek, 
& Morton 2011). The allowed range of Alfven wavelengths 
are bounded on the large scale by the size of the cloud 
and on the small scale by the ion-neutral mean free path 
(Mclvor 1977; Carlberg & Pudritz 1990). Observationally, 
the Alfvenic character of molecular clouds is ill-defined. Re- 
cent compilations of Zeeman measurements of the HI 21cm 



line in absorption against continuum sources by Heiles & 
Troland (2005) and OH emission toward molecular cloud 
cores by Troland & Crutcher (2008) found Ma = ctv/va ~1- 
2. However, optical and infrared polarization studies typi- 
cally measure a uniform field geometry in cloud envelopes 
(Moneti et al. 1984; Heyer et al 1987; Goodman et al. 1990), 
which indicates a field resistant to entanglement by strong 
turbulent eddies (Chandrasekhar & Fermi 1953) and sub- 
Alfvenic turbulence. 

An indirect measure of the influence of the interstellar 
magnetic fleld upon the cold, neutral gas component is the 
degree of velocity anisotropy aligned with the magnetic field 
that can result from MHD turbulence. In this context, ve- 
locity anisotropy refers to the difference in directional power 
spectra or structure functions of the velocity field along 
two orthogonal directions (Goldreich & Sridhar 1995). Such 
magnetically aligned velocity anisotropy can be induced by 
the interactions between three or more Alfven wave pack- 
ets on localized scales corresponding to fluctuations of the 
mean fleld component, SB, as proposed by Goldreich & Srid- 
har (1995). On larger scales, a strong mean field threading 
an interstellar cloud may also regulate the flow of material 
(Cho & Lazarian 2003). Computational simulations demon- 
strate that magnetically aligned velocity anisotropy occurs 
only when the fleld is strong (c^/«^=0.01) and the Aflvenic 
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Mach number is less than unity, {Ma < 1) (Cho & Lazar- 
ian 2003; Vestuto, Ostriker, & Stone 2003). More recently, 
Esquivel & Lazarian (2011) used velocity centroids derived 
from numerical MHD simulations to further confirm a rela- 
tionship between velocity anisotropy and gas magnetization 
(B, Ma). They also inferred a dependence of the measured 
anisotropy with the sonic Mach number for highly magne- 
tized gas. 

Using the simulations of Ostriker, Stone, & Gammie 
(2001), Heyer etal (2008) developed and calibrated an anal- 
ysis technique (sec §2) that derives pscudo velocity structure 
functions along two orthogonal axes from imaging data of 
spectral line emission. Transforming the set of velocity and 
density fields generated by the MHD simulations to model 
observations of ^^CO and "CO J=l-0 emission and ap- 
plying their analysis to the model data set, they verified 
that observational signatures of magnetically aligned veloc- 
ity anisotropy are limited to the regime at which the field is 
strong and the motions are sub-Alfvenic. Heyer etal (2008) 
then applied this method to a subfield of imaging observa- 
tions of ^^CO and ^''CO J=l-0 line emission from the Tau- 
rus molecular cloud (Goldsmith etal 2008; Narayanan etal 
2008). They identified velocity anisotropy that is aligned 
with the local magnetic field direction. Assuming a mean 
density and gas temperature and comparing the measured 
degree of anisotropy with values collected from model ob- 
servations derived from the computational simulations, they 
derived a magnetic field strength of 14 /uG and a subcritical 
(Egas/Ecr-it < 1) cloud envelope, where Egns is the cen- 
tral gas surface density and T^crit = {-B)/(63G)^''^ is the 
surface density that parametrizes magnetic support of the 
gas against self-gravity with mean field (B) (Mouschovias & 
Spitzer 1976). 

The sub-field within the Taurus cloud analyzed by 
Heyer etal (2008) is not an unbiased target. The field was 
selected for this analysis owing to the presence of striations 
of i^CO emission that are parallel with the local magnetic 
field direction as derived from optical polarization of back- 
ground stars that assumes the alignment of elongated grains 
with the field. Correspondingly, there was a reasonable ex- 
pectation that the magnetic field plays some role in the gas 
distribution and kinematics in this selected region of the 
cloud. Given this biased and limited area, the more general 
influence of the static magnetic field component or magnetic 
wave phenomena on the gas throughout the Taurus cloud 
has not yet been fully assessed. 

To more fully evaluate the role of the interstellar mag- 
netic field and its variation throughout the Taurus molec- 
ular cloud complex, we have extended the study of Heyer 
etal (2008) by analyzing the presence of velocity anisotropy 
within 180 partitioned subfields of the Taurus Survey map. 
Specifically, we are interested in the conditions over which 
the magnetically aligned velocity anisotropy varies within 
the cloud. Such variations can occur as the cloud ionization 
changes from the envelope to the more deeply-embedded fil- 
aments and core regions or via a transition from sub-Alfvenic 
to super- Alfvenic turbulence. In §2, we briefly summarize the 
molecular line data and the basic analysis, both of which are 
respectively described in more detail in earlier publications 
(Narayanan etal 2008; Heyer etal 2008). In §3, we apply this 
analysis to the partitioned Taurus fleld for both ^^CO and 



CO data sets. The results are placed into a broader phys- 
ical context in §4. 



2 THE DATA AND ANALYSIS METHOD 

AH molecular line data used in this study are part of the 

FCR AO Survey of the Taurus Molecular Cloud described by 
Narayanan etal (2008) and Goldsmith etal (2008). The sur- 
vey imaged ^^CO and ^'^CO J=l-0 emission from 96 deg^ 
of the Taurus cloud with the FCRAO 14 meter telescope. 
The frontend consisted of the 32 element focal plane array 
receiver, SEQUOIA that fed a set of autocorrelation back- 
end spectrometers configured with a spectral resolution of 
25 kHz or ~0.06 km s~^ for both observed transitions. The 
Taurus data used in this study have been further processed 
to account for contributions and losses from the antenna 
error beam (Pineda et al. 2010). The result of this deconvo- 
lution is a calibration of antenna temperatures that, when 
smoothed to the same angular resolution, are congruent with 
values from the CfA 1 meter telescope for which such error 
beam contributions are minimal. In addition, the ^^CO and 
^■'CO spectra have been resampled onto the same spectral 
grid with 0.1 km spacing. 

The calculation of pseudo-structure functions along two 
perpendicular axes is described by Heyer etal (2008). A set 
of eigenvectors and eigenprojections are derived along a se- 
(jueuce of position-velocity slices of a spectroscopic data 
cube using Principal Component Analysis along each car- 
dinal direction of the spectroscopic data cube. Such a de- 
composition measures the spatial variation of line proflles 
along the axis as these differences generate variance. From 
the eigenvectors and eigenprojections that account for sig- 
nal variance above values expected from thermal noise, one 
can respectively derive the velocity difference and scale over 
which the difference occurs in order to build a pseudo- 
structure function. We refer to these measures as pseudo- 
structure functions because our analysis does not directly 
sample velocity differences from the true 3 dimensional ve- 
locity field but rather, extracts velocity displacements from 
the differences in the spectral line profile shapes that are 
an integration of density and projected velocities along a 
line of sight through the cloud. The pscudo structure func- 
tions axe derived for the two cardinal directions of the data 
cube. The velocity anisotropy is measured from the differ- 
ences between the power law parameters that characterize 
the pseudo-structure functions on each cixis. To find the an- 
gle of maximum velocity anisotropy, Oma, the data cube is 
rotated through a sequence of angles, from to 180 degrees 
spaced by 10 degrees. For each angle, the measure of veloc- 
ity anisotropy, ^'(&), is calculated. This modulation allows 
one to fit for the pheise, 6m a and amplitude, *o, for which 
the anisotropy is a maximum. The value of 6m a can be com- 
pared to the local magnetic field direction, ^b, as determined 
by optical and infrared polarization of background stars. 



3 RESULTS 

To assess the degree of velocity anisotropy throughout the 
Taurus Molecular Cloud, we have partitioned the field sur- 
veyed by Narayanan etal (2008) into 180 blocks, each with 
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256x256 pixels. This size corresponds to 1.42x1.42 deg^ or 
3.5x3.5 pc^ at a distance of 140 pc. The blocks are spaced by 
128 pixels in both the R.A. and Declination directions. The 
selected size of each partition is a compromise between hav- 
ing a field sufficiently large to derive useful statistics along 
any axis while limiting the area over which significant vari- 
ations of the magnetic field geometry may occur, which can 
dilute any signature to MHD velocity anisotropy (Cho & 
Lazarian 2003). 

For each partition, the ^^CO and ^''CO emission are 
analyzed for velocity anisotropy following the procedure de- 
scribed by Heyer etal (2008). In this study, we exclusively 
use the normalized differences between the power law am- 
plitudes as a measure of velocity anisotropy as these param- 
eters are more reliably determined in the fits to the pseudo- 
structure functions than the power law indices given the lim- 
ited spatial dynamic range afforded by the position- velocity 
images. The results of this analysis are shown in Figure [T] 
and Figure [2] that show the variation of the anisotropy in- 
dex, '^{9) with 9 and the perpendicular pseudo-structure 
functions evaluated at the angle of maximum anisotropy, 
9 = 9 MA, for '^■^CO and ^'^CO respectively. Strong veloc- 
ity anisotropy is identified as a strong modulation of "^(9) 
with 9 and displaced pseudo-structure functions along the 
two perpendicular axes evaluated ai 9 = 9 ma- Partitions for 
which the calculations are not well defined owing to a lim- 
ited spatial dynamic range with line emission are not shown. 
Such conditions are realized in partitions where there is little 
or no detected emission or when the partition is positioned 
on the edge of an emission feature that truncates the avail- 
able scales to analyze. Since the 12 CO emission extends over 
larger scales than ^^CO emission owing to a larger abun- 
dance within regions of low extinction and radiative trap- 
ping, there are more partitions with valid measures of the 
velocity anisotropy. 

The strong velocity anisotropy identified by Heyer etal 
(2008) in the region of emission striations in the northeast 
sector of the image is confirmed within the individual par- 
titions that comprise this area. In comparison with the rest 
of the cloud, this area exhibits the largest anisotropy ampli- 
tudes. Other notable regions with strong velocity anisotropy 
include partitions along the northern and western edges of 
the Taurus cloud. 

For both ^^CO and ^"^00 , velocity anisotropy is only 
detected in regions of low surface brightness emission. To 
link this anisotropy to the interstellar magnetic field that 
threads the cloud, such anisotropy should be approximately 
aligned with the local field direction. To assess the degree of 
alignment, we have examined the local magnetic field direc- 
tion in each partition using the set of infrared and optical 
polarization vectors (Moneti et al. 1984; Heyer et al. 1987; 
Goodman et al. 1990; Whittet et al 1992). For each parti- 
tion block, we calculate a mean polarization angle, < 9b >, 
derived by weighting each angle by the degree of polarization 

< 9b T.9,P,/EP, (1) 

where the sum is over all polarization measurements within 
the subfield and 9i and Pi are the polarization angle and per- 
cent polarization for the i*'' measurement within the parti- 
tion. A minimum of two polarization measurements within 
a subfield is required to derive < 9b >■ Figure [3] shows 
the cumulative distribution of the absolute value of the dis- 



placement between the angle of maximum anisotropy and 
the mean local magnetic field direction, \9ma~ < 9b > \, 
for both ^^CO and ^^CO data sets. Also shown is the cu- 
mulative distribution expected for a uniformly random dis- 
tribution of angle displacements. The angular displacements 
derived from the ^^CO data appear consistent with a sam- 
ple of values drawn from a random distribution. This is not 
surprising since there is little velocity anisotropy recovered 
by the ^'^CO data so the corresponding angle of maximum 
anisotropy is not well defined. The Kolmogorov-Smirnov (K- 
S) Test provides a quantitative measure of the likelihood 
that two distributions are drawn from the same population. 
The Kolmogorov-Smirnov D statistic of the ^^CO angu- 
lar displacement and a random distribution is 0.162 corre- 
sponding to a probability of 36% that the two distributions 
are similar. In contrast, the distribution of angular displace- 
ments derived from the ^^CO data shows an enhancement at 
small values of \9ma— < 9b > \- Comparing these values to 
a random distribution, the Kolmogorov-Smirnov D statistic 
is 0.292, which implies a very low probability (P=3xl0~^ 
) that the ^^CO derived displacements are drawn from a 
random distribution. 

Figure [T] and Figure [2] illustrate that signatures to ve- 
locity anisotropy are exclusively limited to regions of low 
surface brightness ^^CO and ^^CO line emission. The low 
CO surface brightness emission corresponds to regions of 
low gas column density or visual extinction. The faint sig- 
nal for both 12 CO and ^^CO from these lines of sight is 
due to subthermal excitation conditions corresponding to 
gas volume densities much lower than the critical density of 
CO (n <10'^ cm~"^). To more directly establish this connec- 
tion, the variation of the amplitude of velocity anisotropy, 
^'o, with visual extinction, A^, in each partition is derived 
and shown in Figure U The visual extinctions are derived 
from the 2MASS image presented by Pineda etal (2010). 
Typical error bars derived from the fitting errors of 4" (6*) 
for 12 CO and ^''CO are shown in the legend. The extinc- 
tions displayed in Figure |4] are averaged over the 1.4x1.4 
deg2 area of each partition. At this resolution, the dynamic 
range of extinction values is much smaller than the original 
Ay image of Pineda etal (2010). The data show that sig- 
nificant velocity anisotropy is limited to regions with mean 
extinctions less than 2 magnitudes. 



4 DISCUSSION 

Velocity anisotropy aligned with the magnetic field is a ba- 
sic feature of strong MHD turbulence in the interstellar 
medium (Kraichnan 1965; Montgomery & Turner 1981; Hig- 
don 1984; Goldreich & Sridhar 1995; Lithwick & Goldreich 
2002; Cho & Vishniac 2000). For most cases, the anisotropy 
results from the interaction between oppositely propagating 
Alfven wave packets that differentially distribute kinetic en- 
ergy transverse and parallel to the local magnetic field. The 
earliest versions of these descriptions focused on density fluc- 
tuations in a highly ionized, magnetic medium inferred from 
scintillation measures. However, in the regime of strong ion- 
neutral coupling, these processes should also be applicable 
to the cold, mostly neutral molecular ISM as traced by the 
molecular line emission used in this study subject to the 
limitation that the wavelength of any disturbance must be 
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larger than the ion-neutral mean free path (Arons & Max. 
1975). 

The degree of velocity anisotropy predicted by Goldre- 
ich & Sridhar (1995) increases with decreasing scale as eval- 
uated within the local frame of reference over which the 
magnetic field direction is uniform. That is, turbulent eddies 
become more elongated as these interact with the local mag- 
netic field. Due to the effects of the Alfven waves, one would 
expect some degree of curvature and wandering of the field 
lines over various scales corresponding to the wavelengths 
of the disturbances. Ideally, one would evaluate the velocity 
anisotropy within the local frame of reference. Cho & Vish- 
niac (2000) describe two methods to identify the local frame 
for any location within a 3 dimensional field as provided 
by computation siniulatious. Evaluating eddy shapes within 
this frame of reference both Cho & Vishniac (2000) and 
Cho, Lazarian, & Vishniac (2002) confirm the GS95 scaling 
between parallel and perpendicular wave vectors. However, 
such methods arc not applicable to real observations that 
are limited to projected distributions of the magnetic field 
geometry on the plane of the sky. While we have attempted 
to account for projected, large scale curvature of the mag- 
netic field by deriving a local, mean magnetic field direc- 
tion in each partition, these measures are not sensitive to 
deflections of the field geometry along the line of sight that 
must be present owing to the transverse perturbations of the 
Alfven waves themselves. Nor do these account for any pro- 
jected substructure within each 1.4x1.4 deg^ area. Further- 
more, since our method is sensitive to the field component 
projected onto the plane of the sky, any inclination of the 
magnetic field would further reduce the degree of measured 
anisotropy. Yet, despite these geometric limitations, veloc- 
ity anisotropy aligned along these projected, mean magnetic 
field directions is measured. By not accounting for the in- 
clination of the mean magnetic field relative to the plane of 
the sky and the perturbations of the field direction along 
the line of sight, our measures of anisotropy arc necessarily 
lower limits to the true degree of anisotropy. It is not clear 
whether these measures of anisotropy arise from the effects 
described by GS95 or are due to some other processes that 
differentially regulate the cascade of energy with respect to 
the mean field direction within the partition area. 

The comparison of our method to measure velocity 
anisotropy with numerical simulations demonstrates the lim- 
ited regime over which such magnetic asymmetry is present 
(Heyer ctal 2008). This regime is characterized by strong 
magnetic fields ((c/va)^ << 1) and sub-Alfvenic motions, 
CTv/vA < 1. Both conditions lead to magnetic field lines that 
are resistant to being tangled by vortical motions such that 
the Alfven wave induced perturbations to the magnetic field, 
5B, are small with respect to the mean field, Bo. Based on 
this previous calibration effort, we infer that such strong 
magnetic conditions prevail in the regions in the Taurus 
cloud where clear evidence for magnetically aligned veloc- 
ity anisotropy is identified. These regions correspond to the 
molecular envelope of the cloud that is traced by the low sur- 
face brightness ^^CO emission. This component also shows 
evidence for magnetically aligned striations with patterns 
that are consistent with wave phenomena within the cloud. 

Differential behavior of the measured anisotropy is iden- 
tified between ^^CO and ^^CO . The lower opacity ^^CO 
J=l-0 line effectively probes the column density range 4 



< A„ < 10 magnitudes (Goldsmith etal 2008; Pineda etal 
2010). Owing to increased optical depth and chemical de- 
pletion of ^^CO and ^"^CO molecules onto grains within 
the cold, high density core regions, neither molecule reliably 
probes this high volume and column density regime of the 
Taurus cloud. For an optically thin line, low surface bright- 
ness emission generally corresponds to lower column density 
relative to that probed by high surface brightness signal. For 
the optically thick '^^CO line, the faint signal within this low 
surface brightness component is due to subthermal excita- 
tion conditions corresponding to low gas volume densities. 
Given the excitation requirements and the limits imposed 
by opacity and chemical depletion, the range of volume den- 
sities responsible for the observed ^'^CO and ^'^CO J=l-0 
emissions is lO'^ < n < 7500 cm~^. 

We propose that the observed differential measurements 
of velocity anisotropy between the low and high surface 
brightness regions and 12 CO and "CO emission reflect a 
transition of the MHD turbulent flow. In the low volume 
density, low column density regime, the field is strongly 
coupled to the neutral material and the motions are sub- 
Alfvenic. In the high column density regime of the cloud, the 
gas motions arc super- Alfvcnic. The Alfvcnic Mach num- 
ber in the Taurus envelope, Ma,o, is not directly estab- 
lished from our data. Following Heyer etal (2008), we es- 
timate Ma,o ~ 0.5 based on the comparison of the degree of 
anisotropy with the models. Therefore, the Alfvenic Mach 
number must increase by a factor greater than 2 between the 
envelope and dense filaments to account for the differential 
behavior of velocity anisotropy. 

The transition from sub-Alfvenic motions in the cloud 
envelope to super- Alfvenic velocities in the high column den- 
sity interior of the cloud offers insight to the evolution of 
material and the role of the magnetic field. The Alfven ve- 
locity is expressed as va = (B) / •^/^Wjlmipn, where (B) is 
the mean magnetic field strength and n is the gas volume 
density within the volume l'^ . Owing to increasing density 
and the amplification of the magnetic field, va is expected 
to vary between the envelope and the dense, filaments. 

A simple, heuristic scaling relationship is constructed 
that describes the variation of the Alfvenic Mach number. 
Ma with the observed physical properties of the cloud. The 
gas density, n can be expressed n = "S/l, where S is the 
gas mass surface density and / is the size scale of the vol- 
ume, l^. Both E and I are more accessible with observational 
methods than n. The magnetic field, (B), scales with den- 
sity as (B) = Bo{n/no)"' ■ Combining these relationships, 
the Alfven velocity scales as 

VA ~ (S)/ni/2 = VAA^/^oT-^'^ilo/iy-^'^ (2) 

where va,o is the Alfven velocity with conditions Eo at size 
lo- Any geometrical factors related to re-casting the volume 
density to E/Z are folded into va,o- The velocity disper- 
sion characterizing clouds motions within a volume, is 
described by the velocity structure function, 

a, = a^fi{T./T.oY'\l/loy' (3) 

We include the dependence of the velocity structure func- 
tion on the gas surface density within a GMC. Heyer ct al. 
(2009) identified a link between the gas surface density of a 
GMC and av,L/L^^'^, where crv,L is the velocity dispersion 
measured over the cloud size, L. If the power law index of 
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the first order structure function, 7 = 1/2, as determined 
by Heyer & Brunt (2004), then the quantity au,L/L^^'^ pro- 
vides the normalization of the velocity structure function 
of the GMC. Therefore, we necessarily include the factor, 
(E/Eo)'''^ in the expression for cr„. For the remainder of 
this paper, we will assume that 7 = 1/2. The Alfvenic Mach 
number, Ma then scales as 

Ma = o^/vA = a,,o(E/Eo)''''(V^o)'''V^'A (4) 

Ma = Ma,o(E/So)^-'=(V^o)" (5) 

where Ma,o is the Alfvenic Mach number at conditions with 
surface density Eo and size scale, lo- 
in Figure O we show the variation of MA/MAfi with 
size scale and surface density for three physical cases k = 
0, 1/2, 2/3. The value of ft = 2/3 corresponds to a spheri- 
cally contracting cloud with strict flux freezing conditions; 
K = 1/2 is expected from an isothermal, gravitationally con- 
tracting fragment or core with frozen flux and mass to flux 
ratio greater than the critical value, {M/<^)cr ~ 1/(63G)^''^ 
(Mouschovias & Spitzer 1976). For material falling down 
magnetic field lines as in the case of ambipolar diffusion over 
the limited density range probed by '^^CO and ^^CO J=l- 
emission (Fiedler & Mouschovias 1993) or shocks aligned 
with the mean field field direction, k ~ (Hennebelle et 
al. 2008; Heitsch, Stone, & Hartmann 2009). It is notewor- 
thy that for n — 0, Ma/Ma.o has no dependence on size 
scale (horizontal lines) so that Ma/Ma,o increases with any 
increase in surface density as expected between the cloud 
envelope and dense cores or filaments. For the other cases 
(k = 1/2, 2/3), a much higher contrast of mass surface den- 
sity (E/Eo) over a limited change in size scale is required 
to significantly increase Ma/Ma,o- The amplification of the 
field with increasing density constrains va/va,o 5S 1 while 
o'v/<Jv,o decreases due to the energy cascade of turbulence. 

A visual inspection of the Taurus maps provides a 
coarse estimate of the conditions between the cloud enve- 
lope and the dense filaments: l/lo ^ 1/3 and E/Eo ^ 3. For 
these conditions, the Alfvenic Mach number would remain 
constant between the envelope and dense filament regimes 
for K — 1/2 and would even decrease for k — 2/3. These 
observational constraints exclude magnetic conditions for 
which K 1/2. In contrast, the accumulation of material 
along magnetic field lines without significant magnetic field 
amplification (k ~ 0) can readily achieve the necessary fac- 
tor, Ma/Ma,o ^ 2, within the observed conditions of the 
Taurus cloud. 

Near zero values of k over the density ranges corre- 
sponding to our ^^CO and ^"^CO observations are predicted 
for molecular clouds whose evolution is regulated by ambipo- 
lar diffusion or the large scale shocks perpendicular to the 
field that forces material along the mean field direction. The 
latter scenario could occur as two clouds moving in opposite 
directions along the magnetic field collide (Heitsch, Stone, 
& Hartmann 2009). While this description is plausible, the 
requirement that the relative motions be strictly along the 
magnetic field would make this an improbable event. Even 
a small, angular offset between the ordered magnetic field 
direction and that of the converging fiow would readily re- 
sult in a field component parallel to the compressed layer 
that would be further amplified by the local dynamic pres- 
sure of the shock. Alternatively, ambipolar diffusion is an 



inevitable process in the regime of low ionization and an 
initially sub-critical cloud envelope. A recent polarization 
study by Chapman et al. (2011) derived sub-critical condi- 
tions for several of the major subregions of the Taurus cloud 
over scales of ~ 2 pc. We propose that the observed veloc- 
ity and column density structure of the Taurus cloud is a 
direct consequence of ambipolar diffusion as this provides a 
natural explanation for the differential behavior of velocity 
anisotropy between the cloud envelope and the high column 
density filaments. 



5 CONCLUSIONS 

We have analyzed ^^CO and ^"'CO J=l-0 emissions within 
equally partitioned areas of the Taurus molecular cloud for 
velocity anisotropy. Significant anisotropy aligned with the 
local magnetic field direction is limited to the envelope of 
the cloud as traced by low surface brightness ^^CO emis- 
sion. The "CO emission, which effectively probes the dense 
filaments, exhibits little evidence for such anisotropy. This 
differential behavior between the low column density enve- 
lope and the dense interior results from a transition from 
sub- Alfvenic turbulence to super-Alfvenic motions. Such a 
transition is realized if the magnetic field is not amplified 
with increasing density as expected from ambipolar diffu- 
sion or forced loading of field lines by magnetically aligned 
cloud-cloud collisions. 
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Figure 1. (top) The variation of the anisotropy index, ^(d) with rotation angle of the data cube for each sub-field with sufficient signal 

to carry out the analysis using the -'^■^CO J=l-0 data. The solid line shows the fit of the points to the function 't = ^oCos[2((? — QMAi\- 
(bottom) The pseudo velocity structure functions along the cardinal axes (solid circles; x axis; open circles; y axis) evaluated at the angle 
of maximum anisotropy, 6m A- 
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Figure 2. Same as Figure [T] with '^^CO J=l-0 data. 
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Figure 3. The cumulative distribution of the absolute value of the difference between the angle of maximum anisotropy and the local 
magnetic field direction for ^-^CO (histogram), '^^CO (dotted line), and a random distribution of angle differences (solid line). 
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Figure 4. (left) The variation of the amphtude of velocity anisotropy, 'I'o witli the mean visual extinction over each subfield. The solid 
and open circles are derived from ^^CO and ^^CO respectively. Magnetically aligned velocity anisotropy is only evident within the low 
column density envelope regions of the Taurus cloud. 
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Figure 5. Contours of Ma/Ma,o values as a function of size scale, / and surface density, S normalized to initial values /o and So 
for K = (left panel), 1/2 (middle panel), and 2/3 (right panel). The circumscribed symbols approximate the relative size scale and 
surfax;e density conditions between the low column density envelope (E) and the dense filaments (F) of the Taurus cloud. A transition of 
sub-Alfvenic to super-Alfvenic motions {Ma/Ma,o > 2) between the cloud envelope and dense filaments is only possible for k = 0. 



